Project presentation

Laurine Panel, Rasmus Tøffner-Clausen, Mathias Glerup Filtenborg, Alexander Videbæk, Qiju Chen

2025-12-02

Introduction

Dataset

The dataset focuses on patients with prostate cancer : one of the most common cancer types among men
Prostate dataset from this data base : https://hbiostat.org/data/

Aim

Find out which treatment is the most effective
Determine which indicators are significant
Apply what we have learned

Materials and methods - Data preparation

Downloaded into two separate sets:

Metadata

data_dir <- "../data/"
meta_tsv_file <- "01_meta_dat_load.tsv"
if (!file.exists(str_c(data_dir, meta_tsv_file)) ) {
  prostate |>
    select(patno, stage, status, age, wt, pf, hx, sdate) |>
    write_tsv(file = str_c(data_dir, meta_tsv_file))

Treatment data

treatment_tsv_file <- "01_treatment_dat_load.tsv"
if (!file.exists(str_c(data_dir, treatment_tsv_file)) ) {
  prostate |>
  select(-stage, -status, -age, -wt, -pf, -hx, -sdate) |>
    write_tsv(file = str_c(data_dir, treatment_tsv_file))

Then converted the two raw files into a tidy dataset:
- by removing NA values
- creating survival event variables (0 and 1)
- removing ap outlines (exceeding 100 to avoid distorsion)

Materials and methods

For modelling we created two augmentations.

df_long <- df_num |>
  select(c(patient_no,stage, age, weight, follow_up_months,
                        sys_blood_pressure, dia_blood_pressure,
                        serum_hemoglobin, serum_prostatic_acid_phosphatase,
                        combined_index, size_of_primary_tumor, is_alive, performance, medical_history,ekg, bone_metastases, treatment, dosage)) |> 
  pivot_longer(cols = c(stage, age, weight, follow_up_months,
                        sys_blood_pressure, dia_blood_pressure,
                        serum_hemoglobin, serum_prostatic_acid_phosphatase,
                        combined_index, is_alive, performance, medical_history,ekg, bone_metastases, treatment, dosage),
               names_to = "variable",
               values_to = "value")

augment_data <- full_df |>
  mutate(MAP = (sys_blood_pressure+ 2 *dia_blood_pressure) / 3,
         is_alive = case_when(is_alive == "alive" ~ 0,
                              is_alive == "dead" ~ 1),
         treatment_detail = case_when(
           treatment == "placebo" ~ "Placebo",
           treatment == "estrogen" & dosage == 0.2 ~ "0.2mg Estrogen",
           treatment == "estrogen" & dosage == 1.0 ~ "1.0mg Estrogen",
           treatment == "estrogen" & dosage == 5.0 ~ "5.0mg Estrogen",
           TRUE ~ NA_character_),
         treatment_detail = factor(treatment_detail, 
                                   levels = c("Placebo", 
                                              "0.2mg Estrogen", 
                                              "1.0mg Estrogen", 
                                              "5.0mg Estrogen")),)

Results - Preliminary model

Linear Volcano

Linear Tree

Results - Preliminary model

QQ-plot linear

QQ-plot log

Results - Preliminary model

Log-Volcano

Log-Forest

Results - Survival Analysis

Significant risk factor : BM, tumor size, AP and age
Significant protective factor : hemoglobin and 1.0 mg estrogen

Project presentation

Introduction

Dataset

Aim

Materials and methods - Data preparation

Metadata

Treatment data

Materials and methods

Results - Preliminary model

Results - Preliminary model

Results - Preliminary model

Results - Survival Analysis

Results - Survival Analysis

Results - Survival Analysis

Discussion

How we can improve this project ?

Tons more data we’d like to model!

Building our own prediction model based on AP-levels.

Comparison with antigen data!